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ABSTRACT 

Rotational excitation of isotopologues of HCN and HNC by thermal electron- 
impact is studied using the molecular R-matrix method combined with the adiabatic- 
nuclei-rotation (ANR) approximation. Rate coefficients are obtained for electron tem- 
peratures in the range 5—6000 K and for transitions among all levels up to J = 8. 
Hyperfine rates are also derived using the infinite-order-sudden (IOS) scaling method. 
It is shown that the dominant rotational transitions are dipole allowed, that is those 
for which A J = 1. The hyperfine propensity rule A J = AF is found to be stronger 
than in the case of He— HCN collisions. For dipole allowed transitions, electron-impact 
rates are shown to exceed those for excitation of HCN by He atoms by 6 orders of 
magnitude. As a result, the present rates should be included in any detailed popula- 
tion model of isotopologues of HCN and HNC in sources where the electron fraction 
is larger than 10 -6 , for example in interstellar shocks and comets. 



1 INTRODUCTION 



Hydrogen cyanide, HCN, and its isomer HNC are among 
the most abundant organic molecules in space, from star- 
forming regions to circumstellar envelopes and comets. HCN 
and HNC also belong to the small class of molecules detected 
in high-redshift galaxies, along with CO and HCO + (Guelin 
et al. 2007 and references therein). In addition to thermal 
emission from various rotational transitions within different 
vibrational states at (sub)millimeter and far-infrared wave- 
lengths (see, e.g., Cernicharo et al. 1996), a few masering 
lines have been detected toward several stars (e.g. Lucas & 
Cernicharo 1989). Numerous isotopologues have also been 
identified, in particular the deuterated species DCN and 
DNC (e.g Leurini et al. 2006 and references therein). 

One interesting observation concerning rotational spec- 
tra of HCN is that of the 'hyperfine anomalies' (Walmsley, 
Churchwell & Fitzpatrick 1982 and references therein). At 
high resolution it is possible to resolve the hyperfine compo- 
nents arising from the nitrogen ( 14 N) nuclear spin for transi- 
tions arising from low-lying rotational levels. The hyperfine 
lines have been found in a number of cases (e.g. Izumiura, 
Ukita & Tsuji 1995; Park, Kim & Minh 1999; Ahrens et 
al. 2002 and references therein) to be not in thermal equi- 
librium with each other. These anomalies have been shown 
to depend on the degree of thermal overlap, on the opacity 
and on the collisional rates (e.g. Guilloteau & Baudry 1981; 
Gonzalez- Alfonso & Cernicharo 1993; Turner 2001). 

Already thirty years ago it was suggested that when 
the electron fraction exceeds ~10 -5 , electron-impact exci- 



tation of polar molecules may be significant in addition to 
collisions with the most abundant neutrals (Dickinson et 
al. 1977). Such conditions can be found for example in dif- 
fuse interstellar clouds and photon-dominated regions where 
n(e)/n(H) can reach a few 10 -4 . Recently Lovell et al. (2004) 
devoted their study to the effect of electrons in the rotational 
excitation of cometary HCN. They showed that electron col- 
lisions are the dominant mechanism of HCN excitation in the 
comets Hale-Bopp and Hyakutake where the electron frac- 
tion, n(e)/n(H2 0), lies in the range ~ 10~ 5 -1. The authors 
stated that accounting for electron collisions may thus alle- 
viate the need for large HCN-H2O cross sections in models 
that neglect effects due to electrons. Similar conclusions were 
drawn in the case of water (Faure, Gorfinkiel & Tennyson 
2004a and references therein). 

Recently Jimenez-Serra et al. (2006) used rotational 
emissions of SiO, HCO + , HCN and HNC to probe electron 
densities in C-type shocks. In particular they ascribed differ- 
ences in the ambient and precursor components to electron 
density enhancements during the first stages of the C-type 
shock evolution. The dipolar ion HCO + was found to be a 
sensitive tracer of this effect. 

A theoretical determination of rotational cross sec- 
tions and rate coefficients for e— HCN was carried out by 
Sana et al. (1981) for the rotationally inelastic transitions 
J = — > 1, J = 1 — > 2 and J = — + 2. Cross sections 
were computed using the rotational close-coupling method 
and rate coefficients were obtained in the temperature range 
5—100 K. An additional study includes that of Jain & Nor- 
cross (1985), where the adiabatic-nuclei-rotation (ANR) ap- 
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proximation was combined with model potentials. The cross 
sections were computed for the same set of transitions as 
in Saha et al. (1981). In addition to electron scattering, the 
rotational excitation of HCN by He as substitute for H2 
has been studied by Green & Thaddeus (1974), Green (un- 
published data J] and Monteiro & Stutzki (1986) (see also 
references therein) . This latter work also provided hyperfine 
rates for J=0-4 and T=10-30 K. 

In the present study, the molecular R-matrix method 
has been combined with the adiabatic-nuclei-rotation 
(ANR) approximation. We provide e— HCN rotational rates 
for transitions among all levels up to J=8 and for tempera- 
tures in the range 5—6000 K. Hyperfine rates are obtained 
within the infinite order sudden (IOS) scaling method for 
J = — 3. Comparisons to the previous theoretical studies 
and to differential data is also presented. In Section 2, R- 
matrix calculations are described and the procedure used to 
obtain rotational and hyperfine rates is briefly introduced. 
In Section 3, both cross sections and rate coefficients are 
presented and discussed. Conclusions are given in Section 4. 



2 CALCULATIONS 
2.1 R-matrix calculation 

Two of us (Varambhia & Tennyson 2007) have carried 
out R-matrix calculations on HCN and HNC at the close- 
coupling level in C21) symmetry. The calculations were per- 
formed using the UK molecular polyatomic R-matrix pack- 
age (Morgan, Tennyson & Gillan 1998) at the (fixed) equilib- 
rium geometries of HCN and HNC using an R-matrix radius 
10 a . 

HCN and HNC (hence DCN and DNC) were repre- 
sented by the 6-31G Gaussian type orbital (GTO) basis 
set. The target wavefunctions were computed using the 
complete active space configuration interaction (CASCI) 
method. They were subsequently improved using a pseudo 
natural orbitals calculation. The five lowest C2i> electroni- 
cally excited states 1 Ai, 3 Ai, 3 A2, 3 Bi, 3 B2 were employed 
and all possible single and double excitations to virtual or- 
bitals were included. In order to incorporate the double ex- 
citation however, it was necessary to freeze eight electrons 
(the Is and 2s electrons of the C and N atoms). For both 
HCN and HNC the weighting coefficients for the density ma- 
trix averaging procedure were 5.75, 1.5, 1.5, 1.5, 1.5 for J Ai, 
3 Ai, 3 A 2 , 3 Bi, 3 B 2 respectively. This target model yielded 
dipole moments -2.87 D and -2.91 D for HCN and HNC 
respectively, which can be compared to the experimental 
values -2.985 D (Ebenstein & Muenter 1984) and -3.05 D 
(Blackman et al. 1976). 

The scattering wavefunctions were represented by a 
close-coupling expansion including 24 target electronic 
states, done to keep the expected shape resonance as low 
as possible. Calculations were performed on the scattering 

1 The earlier calculations of Green & Thaddeus (1974) (re- 
stricted to the lowest 8 levels of HCN and temperatures of 
5—100 K) were extended in 1993 to obtain rate constants 
among the lowest 30 rotational levels and for temperatures 
of 100—1200 K. These unpublished results are available at 
http : //data. giss .nasa.gov/mcrat es/#hcn 



states 2 Ai, 2 Bi, 2 B 2 , 2 A 2 . The continuum GTOs are those 
of Faure et al. (2002) and are expanded up to g (1=4) partial 
wave. One virtual orbital was allocated to each symmetry 
where target orbitals were available to do so. This scattering 
model also used the separable treatment of continuum and 
virtual orbitals in order to allow for increased polarizabil- 
ity and hence further reduce the position of any resonance 
yielded by the model. Further discussion on the problem of 
convergence of the polarisation interaction in close-coupling 
methods can be found in Varambhia & Tennyson (2007), 
Gil et al. (1994) and Gorfinkiel and Tennyson (2004). The 
resulting fixed-nuclei (body fixed) T-matrices were used to 
obtain the pure rotation excitation cross-section in the col- 
lision energy range 0.01—6.2 eV (see Section 3.2). The posi- 
tions of the shape resonances for these molecules are inside 
the scattering energy ranges. 

In the fixed-geometry approximation, the DCN and 
DNC electronic wavefunctions are identical to those of HCN 
and HNC, respectively. The distinction between the isotopo- 
logues therefore arises, in the present treatment, only from 
the different rotational excitation thresholds (see below). 



2.2 Rotational cross sections 

We consider here an electron scattering from a polar lin- 
ear neutral molecule in the fixed-nuclei (FN) approxima- 
tion (Lane 1980). In this approach, the cross sections are 
expressed as a partial-wave expansion within the ANR ap- 
proximation which assumes that the initial and final target 
rotational states are degenerate. For low partial waves (here 
I <4), the cross section is computed from the FN T-matrices 
obtained via the R— matrix calculations. In the case of dipole 
forbidden transitions (i.e. those with A J > 1), cross sections 
are expected to converge rapidly and can be safely evaluated 
using FN T-matrices only (see Section 3.1). In the case of 
dipole allowed transitions (A J = 1), however, the partial- 
wave expansion does not converge in the FN approximation, 
owing to the very long-range nature of the electron-dipole 
interaction. To circumvent this problem, the standard proce- 
dure is to use the dipolar Born approximation to obtain the 
cross section for the high partial-waves not included in the 
FN T-matrices (Crawford & Dalgarno 1971). In this case, 
the final cross section is calculated as the sum of two con- 
tributions and can be regarded as a 'short-range correction' 
to the (dipolar) Born approximation. 

The known unphysical behaviour of the FN cross- 
sections near rotational thresholds, inherent in the ANR 
approximation, is corrected using a simple kinematic ratio 
(Chang & Temkin 1970) which forces the excitation cross 
sections to zero at threshold. In the case of e— H2, this pro- 
cedure has been shown to be accurate down to a collision 
energy E ~ 2 x AE where AE is the rotational threshold 
(Morrison 1988). Recently, experimental data for the scat- 
tering of cold electrons with water has also confirmed the 
validity of the adiabatic 'threshold-corrected' approximation 
down to very low electron energies (Curfk et al. 2006; Faure, 
Gorfinkiel & Tennyson 2004b). Note that threshold effects 
could only be included rigorously in a full rotational close- 
coupling calculation, which is impractical at the collision 
energies investigated here. In this context, it is worth stress- 
ing that in the case of molecular ions, cross sections are 
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large and finite at threshold with a significant but moderate 
contribution from closed-channels (Faure et al. 2006) . 

Finally, in order to make comparisons with available 
experimental data, differential cross sections were computed 
using the Born-closure method for the scattering amplitude, 
as described in Itikawa (2000). In this approach, the high- 
partial waves contribution (I > 4) due to the quadrupolc 
and induced- dip ole interactions was also included. 



k IOS (JF -> J'F') = 



(2J + 1)(2J' + 1)(2F' + 1) 

2 



2.3 Hyperfine rates 

The hyperfine interaction arises due to the very weak cou- 
pling of nuclear spin to the molecular rotation, which to 
an excellent approximation does not affect the overall colli- 
sion dynamics. In the laboratory, Ahrens et al. (2002) mea- 
sured ground state rotational transitions of HCN using sub- 
Doppler saturation spectroscopy in the THz region. This 
technique enables features such as hyperfine structure to be 
revealed, which within Doppler limits, would remain hid- 
den. Nine consecutive rotational transitions with their as- 
sociated hyperfine structure have been thus partly resolved. 
Additional study includes that of Turner (2001). Here the 
transitions J = I — » and 2 — ► 1 were observed toward 
dark clouds for the deuterated molecules N2D + , DCN and 
DNC, from which molecular constants including the com- 
plex nuclear quadrupole hyperfine splitting were derived. 
More recently the nuclear quadrupole hyperfine structure of 
HNC and DNC was resolved in the laboratory for the first 
time using millimetre-wave absorption spectroscopy (Bech- 
tel, Steeves & Field 2006). New rest frequencies for the 
J — 1 — > 0, J = 2 ^ 1 and J = 3 — > 2 rotational transi- 
tions of the ground vibrational state were determined. It was 
found that the hyperfine structure of HNC is dominated by 
the interaction of the valence shell electrons with the nuclear 
spin of the nitrogen atom. The hyperfine structure of DNC, 
however, was much more complicated due to the additional 
coupling of the deuterium nucleus (1=1). This coupling gave 
rise to a septet in the J = I — > transition. 

As hyperfine structure is resolved in some of the astro- 
nomical spectra, it is crucial to know the rate coefficients 
among these energy levels. In the following, only the spin of 
the 14 N atom is considered as the effects of the spins of H 
and D atoms on the collisional rates is assumed negligible, 
although of course they affect the energy levels. Hyperfine 
states are labeled by J and F, where F is the total angular 
momentum obtained by coupling J to /, the nuclear spin 
(1 = 1 for 14 N). If the hyperfine levels are assumed to be 
degenerate, it is possible to simplify considerably the hyper- 
fine scattering problem (Corey & McCourt 1983). Within 
the ANR or IOS approximation, which both ignore the ro- 
tational spacings, the scattering equations are further sim- 
plified and the rates among rotational and hyperfine levels 
can be simply calculated in terms of 'fundamental' rotational 
rates, k IOS (L -f 0): 
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are Wigner-3j and Wigner-6j symbols, respectively. 

However, as the ANR rotational cross sections are cor- 
rected for threshold effects (see above), Eqs. (1) and (2) 
are only moderately accurate for the actual rates. We have 
therefore implemented the 'scaling' method proposed by 
Neufeld & Green (1994) in which the hyperfine rates are 
obtained as a scaling of the rotational rates: 

'" IOS (JF — > J'F') 



k(JF -> J'F') = 



k 



J> 



-k(J^j'), 



(3) 



k ms (J 

using the actual rates k(L — > 0) for the IOS fundamen- 
tal rates. For quasi-elastic transitions, i.e. JF — > JF' with 
F F' , Eq. (2) is used directly, as recommended by Neufeld 
& Green (1994). Note that fundamental rates for down- 
ward transitions were employed as these were found to pro- 
vide better results in the case of He— HCN for which accu- 
rate close-coupling hyperfine rates are available (Monteiro 
& Stutzki 1986). 

Finally, it should be noted that within the IOS approx- 
imation, the allowed transitions or selection rules are deter- 
mined by the Wigner-6j symbol, as expressed in Eq. (2). 
Interestingly, the same selection rules have been obtained 
by Chu (1976) and Varshalovich & Khersonskii (1977) from 
multipole expansion approaches. Radiative (dipolar) selec- 
tion rules are also determined by the above 6j symbol with 
L = 1. 



3 RESULTS AND DISCUSSION 
3.1 Cross Sections 

The comparison with experimental differential cross sections 
is often the only reliable way for testing electron-molecule 
calculations. In the case of polar targets, integral cross sec- 
tions deduced from experiments are indeed strongly depen- 
dent on the extrapolation procedure used to estimate the 
small-angle scattering which cannot be detected experimen- 
tally (e.g. Faure et al 2004b). 

We have computed the elastic (rotationally summed) 
differential cross sections for electron scattering from HCN 
at the electron energy of 5 eV in order to compare our cal- 
culations with the experimental results of Srivastava et al. 
(1978). As shown in Fig. 1, our calculations reproduce the 
experimental data only qualitatively. Quantitatively, our re- 
sults are larger by roughly a factor of 2 at angles 9 < 100°. 
It can be noticed that the dipolar transition 0— >1 dominates 
over all other transitions at both forward and backward an- 
gles. Sharp dips are also observed in the 0— >1 transition at 
about 20, 70 and 120°. We note that these dips, which re- 
flect rainbow scattering (e.g. Korsch & Ernesti 1992), are 
not completely suppressed by the contribution of the other 
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Scattering angle (deg) 



Figure 1. Differential cross sections for vibrationally elastic scat- 
tering (rotationally summed) of electrons by HCN at 5 eV. Exper- 
imental data (full circles) are from Srivastava et al. (1978). The 
present calculation is given by the thick solid line. Other lines 
denote partial state-to-state differential cross sections. 




Electron energy (eV) 

Figure 2. Integral cross sections for rotationally inelastic scatter- 
ing of electrons by HCN as a function of electron energy. Lozenges 
and stars denotes the results of Saha et al. (1981) and Jain and 
Norcrosss (1985), respectively. The present calculations are given 
by the solid thick lines. Dotted lines are extrapolations (see text). 



transitions in the total curve. More experimental study is 
clearly necessary to interpret the difference between theory 
and experiment, as already suggested sometime ago by Jain 
& Norcross (1985). 

Note that the Born completion was found to be crucial 
in the case of the dipole interaction, as expected for very 
polar species. Dipole allowed transitions were thus found to 
be largely dominated by the high-partial waves (I > 4) with 
only a small contribution from the FN T-matrices. In con- 
trast, the Born completion for the quadrupole and induced- 
dipole interactions was found to be negligible: high-partial 
waves were found to increase cross sections for A J = 0, 2 
by less than 2%. This clearly shows that cross sections for 
dipole forbidden transitions converge rapidly with respect 
to the partial wave expansion, as observed in the case of 
molecular ions (e.g. Faure & Tennyson 2001). As a result, 
only the Born-closure approximation for the dipole inter- 
action was implemented in the computation of the integral 
rotational cross sections. 

In Fig. 2, we compare our integral cross sections with 
those of Jain & Norcross (1985) and Saha et al. (1981). Ex- 
cellent agreement is observed down to 0.01 eV. At lower en- 
ergy, we adopted the extrapolation formula of Rabadan et al. 
(1998) (Eq. 1 of their paper) which was calibrated using the 
rotational close-coupling results of Saha et al. (1981). This 
procedure obviously introduces uncertainties in the rate cal- 
culation at temperatures below about 100 K, in addition to 
closed-channels effects which are completely ignored in the 
present treatment. Hence, rate coefficients below 100 K are 
expected to be less accurate than at higher temperatures. 

A large shape resonance can be observed in Fig. 2 at 
about 2.8 eV. It is particularly strong for transitions with 
small cross sections, i.e. those with AJ > 2. This resonance 
(of 2 II symmetry) has also been observed experimentally, 
as discussed in Varambhia & Tennyson (2007). For HNC, 
it appears at a slightly lower position of 2.5 eV. It is to be 



noted that the cross sections for AJ = 1 are much larger 
than for the other transitions due to the dominance of the 
long-range dipole interaction for very polar molecules. 

We finally note that, to our knowledge, there is no data 
(theoretical or experimental) on the rotational excitation of 
DCN, HNC and DNC available in literature. 

3.2 Rotational rates 

As HCN and HNC (hence DCN and DNC) have their first 
electronic excitation threshold at 6.63 eV and 6.20 eV, re- 
spectively, the rotational excitation cross sections were com- 
puted in the range 0.01—6.2 eV and extrapolated down to 
rotational thresholds (see Section 3.1). Assuming that the 
electron velocity distribution is Maxwellian, rate coefficients 
were obtained for temperatures in the range 5—6000 K and 
for transitions among all levels up to J=8. De-excitation 
rates were calculated using the detailed balance relation. 
For later use in modelling, the temperature dependence of 
the downward transition rates fc(T), in units of cm 3 s _1 , was 
fitted to the analytic form (Faure et al. 2004a): 

JV 

log 10 k(T) = a ^ ( 4 ) 

i=0 

where x = 1/T 1//6 , N = 4 and T is restricted to the range 
5—2000 K in order to achieve a fitting accuracy of a few per- 
cent. The coefficients {a^}, in units such that k is in cm 3 s _1 , 
obtained from the fitting procedure are given in tables 1, 2, 
3 and 4 for HCN, DCN, HNC and DNC, respectively. Fit- 
ting coefficients for temperatures above 2000 K and J > 8 
can be obtained upon request to the authors. 

Excitation rates are presented in Fig. 3 for HCN and 
DCN. The small differences with the results of Saha et al. 
(1981) reflect the sensitivity of low temperature rates to the 
near-threshold cross sections. It can also be noticed that 
the smaller rotational thresholds of DCN (with respect to 
HCN) lead to a significant but moderate increase of rates 
below 100 K. As a result, rates for other isotopologues such 
as H 13 CN and H 15 NC are expected to be very similar to 
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Table 1. Coefficients a,; of the polynomial fit, equation (4), to the rotational de-excitation rate coefficients of HCN. These coefficients 
are only valid in the temperature range 5—2000 K. _E up (K) are the upper level energies. 
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-Bup(K) 




ai 


"2 


<J3 


(14 


(1-0) 


4.3 


-7.543 


10.018 


-17.974 


14.328 


-4.169 


(2-0) 


12.8 


-10.923 


22.212 


-52.611 


60.265 


-26.660 


(2-1) 


12.8 


-7.552 


10.772 


-21.515 


19.858 


-7.067 


(3-0) 


25.5 


-13.297 


30.136 


-71.602 


79.988 


-34.340 


(3-1) 


25.5 


-10.830 


22.286 


-52.461 


59.201 


-25.830 


(3-2) 


25.5 


-7.589 


11.377 


-24.257 


24.067 


-9.239 


(4-0) 


42.5 


-15.647 


35.452 


-84.761 


94.380 


-40.338 


(4-1) 


42.5 


-13.222 


30.599 


-73.095 


81.660 


-35.036 


(4-2) 


42.5 


-10.827 


22.664 


-53.590 


60.351 


-26.259 


(4-3) 


42.5 


-7.633 


11.910 


-26.550 


27.495 


-10.976 


(5-0) 


63.8 


-17.755 


36.842 


-88.707 


99.539 


-42.965 


(5-1) 


63.8 


-15.718 


37.455 


-91.762 


104.331 


-45.407 


(5-2) 


63.8 


-13.277 


31.655 


-76.702 


86.569 


-37.456 


(5-3) 


63.8 


-10.815 


22.782 


-53.882 


60.382 


-26.134 


(5-4) 


63.8 


-7.677 


12.402 


-28.613 


30.532 


-12.492 


(6-0) 


89.3 


-17.947 


22.002 


-49.290 


54.613 


-24.371 


(6-1) 


89.3 


-18.379 


44.219 


-114.389 


136.990 


-62.411 


(6-2) 


89.3 


-16.049 


41.262 


-105.105 


123.812 


-55.550 


(6-3) 


89.3 


-13.426 


33.410 


-82.796 


95.273 


-41.916 


(6-4) 


89.3 


-10.888 


23.601 


-56.532 


63.817 


-27.757 


(6-5) 


89.3 


-7.737 


13.016 


-31.039 


34.097 


-14.297 


(7-0) 


119.1 


-17.477 


2.037 


2.640 


-10.257 


8.971 


(7-1) 


119.1 


-19.755 


39.818 


-106.846 


130.991 


-59.153 


(7-2) 


119.1 


-20.167 


61.606 


-172.992 


219.967 


-103.973 


(7-3) 


119.1 


-17.003 


50.704 


-137.567 


170.812 


-79.770 


(7-4) 


119.1 


-13.717 


36.432 


-93.345 


110.640 


-49.901 


(7-5) 


119.1 


-10.947 


24.277 


-58.843 


66.970 


-29.308 


(7-6) 


119.1 


-7.791 


13.571 


-33.249 


37.364 


-15.963 


(8-0) 


153.1 


-9.087 


-41.736 


42.984 


14.247 


-28.504 


(8-1) 


153.1 


-14.491 


-23.527 


84.511 


-119.423 


58.957 


(8-2) 


153.1 


-16.613 


11.279 


-11.836 


-2.489 


6.619 


(8-3) 


153.1 


-17.542 


36.928 


-88.046 


95.525 


-39.360 


(8-4) 


153.1 


-18.534 


64.949 


-184.148 


234.228 


-109.912 


(8-5) 


153.1 


-14.334 


42.448 


-114.006 


140.539 


-65.328 


(8-6) 


153.1 


-11.068 


25.526 


-63.164 


73.166 


-32.507 


(8-7) 


153.1 


-7.823 


13.903 


-34.635 


39.376 


-16.952 



those of HCN and HNC, respectively, since the rotational 
thresholds are only slightly different. 

In Fig. 4, electron-impact rate coefficients for HCN and 
HNC are presented along with rates for excitation of HCN 
by He atoms. It is firstly observed that HCN and HNC have 
similar rates, with typical differences less than 50 % for the 
plotted transitions. Note that the dipolar — > 1 curves are 
almost superposed in the whole temperature range. This re- 
flects the dominance of the dipole interaction. Differences 
larger than a factor of 2 were however found for the smallest 
rates, i.e. those corresponding to transitions with AJ > 2. 
Second, rate coefficients for excitation by He atoms are 
much smaller than electron-impact rates. In particular the 
He— HCN rate for the transition — > 1 is about 6 orders of 
magnitude smaller than the e— HCN rate. Furthermore, the 
propensity rules are very different, with even AJ strongly 
favored over odd AJ in the case of He. This propensity re- 
flects an interference effect related to the even anisotropy 
of the He— HCN PES (e.g. McCurdy & Miller 1977). As a 
result, the difference between e— HCN and He— HCN rates 




Temperature (K) 

Figure 3. Rate coefficients for rotational excitation of HCN and 
DCN by electron-impact as a function of temperature. Lozenges 
denote the HCN calculations of Saha et al. (1981) while the solid 
and dashed lines give the present results for HCN and DCN, re- 
spectively. 
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Table 2. Coefficients a,; of the polynomial fit, equation (4), to the rotational de-excitation rate coefficients of DCN. These coefficients 
are only valid in the temperature range 5—2000 K. _E up (K) are the upper level energies. 



ransition 


£ up (K) 


do 


ai 


0,2 


0,3 


0,4 


(1 - 


0) 


3.5 


-7.513 


9.737 


-16.754 


12.437 


-3.174 


(2- 


0) 


10.4 


-10.931 


22.292 


-52.959 


60.944 


-27.038 


(2- 


1) 


10.4 


-7.527 


10.532 


-20.385 


18.095 


-6.144 


(3- 


o) 


20.9 


-13.275 


29.938 


-71.022 


79.434 


-34.150 


(3- 


1) 


20.9 


-10.799 


22.009 


-51.638 


58.299 


-25.449 


(3- 


2) 


20.9 


-7.561 


11.088 


-22.909 


21.979 


-8.157 


(4- 


0) 


34.8 


-15.559 


34.612 


-81.872 


90.381 


-38.336 


(4- 


1) 


34.8 


-13.189 


30.289 


-72.067 


80.410 


-34.476 


(4- 


2) 


34.8 


-10.794 


22.357 


-52.623 


59.213 


-25.753 


(4- 


3) 


34.8 


-7.592 


11.502 


-24.778 


24.819 


-9.611 


(5- 


0) 


52.1 


-17.594 


35.278 


-83.187 


91.562 


-38.855 


(5- 


1) 


52.1 


-15.554 


35.865 


-86.214 


96.384 


-41.332 


(5- 


2) 


52.1 


-13.199 


30.907 


-74.137 


83.071 


-35.733 


(5- 


3) 


52.1 


-10.807 


22.701 


-53.656 


60.318 


-26.194 


(5- 


4) 


52.1 


-7.643 


12.076 


-27.155 


28.395 


-11.451 


(6- 


0) 


73.0 


-17.202 


14.785 


-24.098 


17.750 


-5.212 


(6- 


1) 


73.0 


-17.737 


38.047 


-92.971 


105.803 


-46.244 


(6- 


2) 


73.0 


-15.705 


37.912 


-93.299 


106.476 


-46.487 


(6- 


3) 


73.0 


-13.236 


31.583 


-76.490 


86.276 


-37.307 


(6- 


4) 


73.0 


-10.843 


23.178 


-55.163 


62.141 


-27.005 


(6- 


5) 


73.0 


-7.677 


12.439 


-28.711 


30.698 


-12.601 


(7- 


0) 


97.3 


-16.289 


-7.875 


30.337 


-38.294 


15.702 


(7- 


1) 


97.3 


-18.107 


24.976 


-59.631 


69.686 


-32.154 


(7- 


2) 


97.3 


-18.405 


45.121 


-117.556 


141.673 


-64.869 


(7- 


3) 


97.3 


-16.011 


41.257 


-105.090 


123.775 


-55.517 


(7- 


4) 


97.3 


-13.372 


33.096 


-81.743 


93.795 


-41.171 


(7- 


5) 


97.3 


-10.838 


23.256 


-55.432 


62.368 


-27.064 


(7- 


6) 


97.3 


-7.728 


12.953 


-30.725 


33.628 


-14.069 


(8- 


0) 


125.1 


-9.469 


-41.064 


53.519 


-24.563 


5.360 


(8- 


1) 


125.1 


-17.037 


-1.195 


15.814 


-33.284 


23.214 


(8- 


2) 


125.1 


-19.699 


39.962 


-107.384 


131.888 


-59.701 


(8- 


3) 


125.1 


-19.905 


59.603 


-166.518 


211.269 


-99.909 


(8- 


4) 


125.1 


-16.738 


48.424 


-129.770 


159.575 


-74.013 


(8- 


5) 


125.1 


-13.600 


35.446 


-89.911 


105.637 


-47.303 


(8- 


6) 


125.1 


-10.913 


24.049 


-58.131 


66.107 


-28.940 


(8- 


7) 


125.1 


-7.765 


13.337 


-32.316 


36.014 


-15.296 



is much larger for dipole allowed transitions than for other 
transitions. For instance, the He— HCN rate coefficient for 
the transition — > 4 at 1000 K is only a factor of 4 smaller 
than the electron-impact rate. 

These large differences between He— HCN and e— HCN 
rotational rates, both in terms of magnitude and propensity 
rules, suggest that the modelling of HCN excitation in astro- 
physical environments must be very sensitive to the relative 
abundance of electrons and neutrals. 



3.3 Hyperfine rates 

Hyperfine rate coefficients were computed from the pure ro- 
tational rates using the formalism proposed by Neufeld & 
Green (1994), as explained in Section. 2.3. Sample results 
are presented in Tables 5, 6, 7 and 8 for HCN, DCN, HNC 
and DNC, respectively. Additional hyperfine data at differ- 
ent temperatures and for higher rotational levels can be ob- 
tained upon request to the authors. The relative propen- 
sities predicted by our calculations at 20 K for a few hy- 
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Temperature (K) 

Figure 4. Rate coefficients for rotational excitation of HCN and 
HNC as a function of temperature. The black and red lines denote 
the present results for electron-impact excitation of HCN and 
HNC, respectively. The green line gives the results of Green & 
Thaddeus (1974) and Green (unpublished data) for the rotational 
excitation of HCN by He atoms. 
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Table 3. Coefficients at of the polynomial fit, equation (4), to the rotational de-excitation rate coefficients of HNC. These coefficients 
are only valid in the temperature range 5—2000 K. F up (K) are the upper level energies. 
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01 A 
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oo. 5yo 
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1 1 O QOQ 
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XC\ 7QO 

— 5U. ( oy 


la a \ 
(6 - 4) 


91.4 


—9.107 


10.635 


—23.005 


26.473 


— 12.585 
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(b - 5) 


91.4 


—7.668 


12.643 
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— oU.oUz 


34.339 


— 14.756 
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IZl.o 


— y. f zo 


xq aii 
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121.8 


-12.497 


-22.208 


88.515 


-132.742 


68.301 


(7-2) 


121.8 


-17.644 


39.253 


-96.958 


104.695 


-39.340 


(7-3) 


121.8 


-18.318 


64.133 


-183.445 


234.541 


-110.659 


(7-4) 


121.8 


-13.698 


39.371 


-109.796 


138.629 


-65.302 


(7-5) 


121.8 


-9.814 


17.433 


-45.443 


55.767 


-26.026 


(7-6) 


121.8 


-7.686 


12.855 


-31.518 


35.688 


-15.428 


(8-0) 


156.6 


-12.532 


-8.671 


-48.477 


123.185 


-76.018 


(8-1) 


156.6 


-9.784 


-56.529 


173.568 


-225.735 


106.060 


(8-2) 


156.6 


-12.407 


-22.181 


87.726 


-131.074 


67.259 


(8-3) 


156.6 


-14.860 


14.808 


-20.242 


5.667 


3.941 


(8-4) 


156.6 


-17.161 


50.763 


-126.624 


131.935 


-45.525 


(8-5) 


156.6 


-14.757 


49.398 


-143.428 


186.181 


-89.164 


(8-6) 


156.6 


-10.025 


19.503 


-52.460 


65.752 


-31.138 


(8-7) 


156.6 


-7.546 


2.952 


7.712 


-20.936 


12.334 



perfine transitions induced by electron-impact are given in 
Table 9 and are compared to He— HCN results. Firstly, it 
may be seen that the present implementation of the method 
of Neufeld & Green (1994), i.e. based on downward funda- 
mental rates, gives a reasonable agreement with the accu- 
rate close-coupling results of Monteiro & Stutzki (1986) for 
He— HCN. Errors on rate values were thus found to range be- 
tween typically 10 and 50 %, suggesting that the IOS scaling 
of rotational rates provides a rather good description of hy- 
perfine rates. Second, it is found that the well known propen- 
sity rule AJ = AF is stronger for electron— HCN than for 
He— HCN. This propensity rule was derived and confirmed 
experimentally by Alexander & Dagdigian (1985) for the 
case of atom-molecule collisions. As discussed by these au- 
thors, the AJ = AF propensity is independent of dynam- 
ics and follows properties of Wigner-6J symbols (see Sec- 
tion 2.3). It is thus expected to become increasingly strong 
as J and J' increase, as indeed observed in Tables 5—8. 
The fact that it is stronger for electrons than for He atoms 
is simply due to the smallness of electron-impact rates for 
transitions with AJ > 2. 

It should be noted that in the case of He— N2H 4 " colli- 



sions, Daniel et al. (2005) have found that the AJ = AF 
propensity does not always apply at low energy, owing to 
the presence of closed-channel resonances. The importance 
of such effects for electron-molecule collisions is however un- 
clear and would deserve further studies. 

It can be noticed in Tables 5—8 that in addition to the 
AJ = AF propensity rule, there is a rigorous selection rule 
at the IOS level: the (J,F = J) -> (J' = 1,F' = 0) transi- 
tions are strictly forbidden, following properties of Wigner- 
6j symbols. Monteiro & Stutzki (1986) have shown that for 
He— HCN this selection rule becomes only a propensity rule 
at the close-coupling level. 

Finally we note that the influence of collisional processes on 
HCN hyperfine anomalies is still a matter of debate in the 
literature (e.g. Gonzales- Alfonso & Cernicharo 1993; Turner 
2001; Daniel, Cernicharo & Dubernet 2006). Future radia- 
tive transfer studies including the present hyperfine rates 
will thus help to further assess the role of collision-induced 
transitions between HCN hyperfine levels. In this context, 
it is important to remember that H2— HCN collisional rates 
are currently obtained by scaling He— HCN rates by the re- 
duced mass ratio. While for para-H2(J=0) this assumption 
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Table 4. Coefficients a, of the polynomial fit, equation (4), to the rotational de-excitation rate coefficients of DNC. These coefficients 
are only valid in the temperature range 5—2000 K. _B up (K) are the upper level energies. 
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is probably qualitatively correct, rates for fi2(J 7^ 0) might 
be significantly different, as found for example in the case of 
H2— H2O (Faure et al. 2007 and references therein). Further- 
more, it is worth stressing that rates for excitation of HCN 
by H atoms and H2O molecules are completely unknown. 



4 CONCLUSIONS 

We have calculated electron-impact rotational and hyper- 
fine (de) excitation rates for isotopologues of HCN and HNC. 
These calculations are based on the molecular R-matrix 
method combined with the (threshold corrected) ANR ap- 
proximation. Our results show that such collisions are es- 
sentially dominated by dipolar transitions owing to the 
large dipole moments of the considered species. However, 
short-range effects are important and were included via R- 
matrix wavefunctions as corrections to the Born approxi- 
mation. Dipole forbidden transitions have thus appreciable 
rates which cannot be neglected in any detailed population 
models of isotopologues of HCN and HNC. fn particular, 
we have shown that electron-impact rates are crucial for 



modelling environments where the electron fraction is larger 
than 10 -6 . In this context, we note the recent suggestion 
by Jimenez-Serra et al. (2006) to use rotational emissions 
of HCO + , HCN and HNC to probe electron densities in C- 
type shocks. The present rates should help to investigate in 
more detail the electron density enhancements expected dur- 
ing the first stages of a C-type shock evolution. Finally, our 
work might help to understand and to model the observed 
variable HCN/HNC ratio in comets (Biver et al. 2006 and 
references therein). 
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Table 5. Hyperfine de-excitation rates in cm 3 s 1 for HCN. Where necessary powers of ten are given in parentheses. 
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Table 6. Hyperfine de-excitation rates in cm 3 s 1 for DCN. Where necessary powers of ten are given in parentheses. 
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Table 7. Hyperfine de-excitation rates cm 3 s 1 for HNC. Where necessary powers of ten are given in parentheses. 
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Table 8. Hypcrfine de-excitation rates cm 3 s 1 for DNC. Where necessary powers of ten are given in parentheses. 
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Table 9. Relative propensities for hyperfine collisional rates as a percentage of the total rates for T = 20 K. Present=rates calculated 
here; He (CC)=closc-coupling He-HCN results of Monteiro & Stutzki (1986); He (IOS)=IOS scaling for He-HCN. 
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